Color Neutral Ground State of 2SC Quark Matter 
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We construct a new color neutral ground state of two-flavor color superconducting quark matter. 
It is shown that, in contrast with the conventionally considered ground state with diquark pairing in 
only one color direction, this new state is stable against arbitrary diquark fluctuations. In addition, 
the thermodynamical potential is found to be lower for this new state than for the conventional one. 
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Recent investigations on the QCD phase diagram have 
discovered a rich diversity of color superconducting quark 
matter phases at low temperatures and intermediate den- 
sities ijy. Particularly interesting are possible imphca- 
tions of these results for the physics of compact stars 
as well as heavy ion collision experiments Q address- 
ing the domain of densities and temperatures where the 
strange quarks are still heavy and confined 0,0] . In both 
systems, the constraint of global color neutrahty has to 
be fulfilled. In addition, the constraint of global electric 
neutrality has to be satisfied when macroscopic objects 
like compact stars are considered and flavor changing pro- 
cesses have enough time to adjust u and d quark chemical 
potentials according to /3-equilibrium. 

In view of the nonperturbative character of QCD, the 
theoretical treatment of hadronic matter at the vicinity of 
the phase transitions for low temperatures and finite den- 
sities is a problem of highest complexity, where rigorous 
theoretical approaches are not yet available and lattice 
QCD simulations are up to now not applicable. There- 
fore one has to rely on effective field-theoretical models of 
interacting quark matter, which are built taking into ac- 
count the symmetry requirements of the QCD lagrangian 
and offer the possibility of dealing with the yet simplified 
interactions in a systematic way. Chiral quark models 
of QCD that adopt a current-current form of the inter- 
action with mesonic and diquark components have been 
particularly useful, since the theories can be bosonized in 
a straightforward way. In the meson-diquark representa- 
tion of these models an effective quantum hadrodynamics 
can be derived 6] , but it has not yet been used to explore 
the QCD phase diagram. A preparatory step for this 
formidable task is the investigation of the mean field ap- 
proximation (MFA), where important progress has been 
recently made within a nonlocal, covariant formuation 
that has been extended even to the study of color super- 
conductivity in the QCD phase diagram J^. 

The two-flavor color superconductivity (2SC) phase 
of quark matter has been first considered in instanton- 
motivated QCD models in , where the question of color 



neutrality was disregarded. However, it was soon realized 
that the 2SC phase in which color symmetry is broken by 
the orientation of the diquark field in one of the color di- 
rections (2SC-fe) entails a mismatch in the quark densities 
of paired and unpaired colors provided that color chemi- 
cal potentials are all equal to each other, /ir = Ms = A'b- 
Therefore, in the 2SC-& state color neutrality requires 
color charge chemical potentials to be readjusted so that 
1^8 = ■^{iJ-r + A*g ~ 2/if,) acquires a nonvanishing value 
while /i3 = f (/^r — l^g) remains zero due to the degener- 
acy of the red and green colors. While this adjustment 
of /ig 7^ has long been considered a proper solution of 
the color neutrality constraint 9], a recent investigation 
of fluctuations around the mean field oriented in the blue 
(0,0,1) direction has revealed the instability of this state 
once color neutrality is required ■ In the present paper 
we investigate the entire space of mean-field orientations 
in order to look for color neutral states which are sta- 
ble against fiuctuations. We find that these correspond 
to color neutral symmetric states (2SC-s) for which the 
condensates are equal in modulus in all the three direc- 
tions of the color space. 

As in Ref. ^3|> consider the simplest version of 
the fiavor SU(2) Nambu-Jona-Lasinio model |TlllTllT3l|. 
extended so as to include the quark-quark interaction 
sector and finite chemical potentials 
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Here rh is the diagonal current mass matrix for light 
quarks, Gs and Gd are coupling constants in color sin- 
glet and anti-triplet channels respectively, and V'ia stands 
for quark fields with fiavor index i = u,d and color in- 
dex a = r,g,b (charge conjugated fields are given by 
i'ia = '''I'^l^i'Ta)- The PauH matrices f = {ti,t2,t^) act 
in flavor space, while and e'^ and e"'^'' are totally anti- 
symmetric tensors in flavor and color spaces, respectively. 
In the present letter we are mainly concerned with the ef- 
fect of color neutrality on the ground state of a two-flavor 
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color superconductor, therefore we will restrict here the 
discussion to the flavor symmetric case and consider the 
extension to electrically neutral matter elsewhere. Re- 
garding quark chemical potentials, the elements of the 
matrix 

/I = diag(/ir,Mff;Mfc>Mr,Ms:A'b) (2) 
can be written as 

^J■r = + Ms/S + /is/S , 

^J■b = A*s/3 - 2^8/3 , (3) 

where is the baryon chemical potential, while /ig and 
fi3 are introduced to ensure color charge neutrality. Our 
aim is to discuss the color superconducting phase of the 
model in the mean-field approximation, which in general 
is characterized by nonvanishing diquark condensates 

A„ =-2Gz5(V;[^z7V^e''''"V',7>> o^^r,g,b. (4) 

For standard values of the diquark coupling. Go = jGs, 
there is no simultaneous chiral symmetry breaking in this 
phase. In addition, since for light quarks the current 
quark masses are significantly smaller than the typical 
values of hb and Aq,, for the purpose of the present study 



we can safely neglect both these small masses and the 
corresponding mesonic mean fields. Within this limit, 
we proceed to calculate the thermodynamical potential 
per unit volume at zero temperature. For convenience we 
perform our calculations in Euclidean space, where the 
thermodynamical potential in MFA is given by 

and the constraint of global color neutrality has to be 
obeyed, i.e. color charge densities should vanish 

Qa = -dQ/dfia^O , a = 3, 8. (6) 

Here we have defined — X]a=r g b ^q' ^^^^ ^^"^ space 
integral is regulated as usual by introducing a sharp 
three-momentum cutoff A. The inverse fermion propa- 
gator M is a 48 X 48 matrix in Nambu-Gorkov, Dirac, 
color and flavor spaces, and can be conveniently written 
as 

with 
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where A^ = 175 A^ and (Gg^) ^ = -[(p4TiMa)74+P-7] 
are 4x4 matrices in Dirac space. After some algebra, it 
is seen that the determinant can be cast into the form 



where 



with 



det M={S+ S- )^ 



(9) 



G±|^Re(C± Cf) A^ + cycl.perm.{r56} (10) 
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With this general expression at hand we can investigate 
the problem of finding the most favored state under the 
constraint of the color neutrality, i.e. Q3 = Qg = 0. 



Since the expression in Ea. (|10|l is totally symmetric 
under cyclic permutations of the three color indices, a cu- 
bic symmetry is expected in the three-dimensional space 
spanned by the three color directions along which the 
magnitudes of the colored diquark condensates are the 
coordinates. In the standard Cartesian representation, 
the condensate vector in color space is given by 



A = Ar 



Co + 



Af, 66 . 



(11) 



Due to the above mentioned symmetry, in what follows 
we will restrict our discussion to the sector defined by 
nonnegative values of the Cartesian coordinates Aq,. In 
this representation, the 2SC-& state is given by A2SC-6 = 
(0, 0, Afc), whereas the color symmetric state is A2SC-S — 
(As, As, As), i.e. a vector pointing from the center to one 
of the edges of a cube in color space. For the 2SC-S state 
color neutrality is achieved with /ig = = 0, while 
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everywhere else color symmetry is broken, entailing that 
either ^g, fj,^ or both have to be different from zero in 
order to fulfill the color neutrality constraint. 

The observation made in 10] stating that a 2SC-6 state 
defines a saddle point of the thermodynamical potential 
^ in the order parameter space, being a minimmn in the 
blue direction but a maximum in the red and green ones, 
leads to the important fact that the 2SC-5 state widely 
considered in the literature is not the true ground state 
of quark matter in the 2SC phase. Now the problem is 
to find the true 2SC ground state, which should be ther- 
modynamically more favorable not only by a lower en- 
ergy but also by its stability with respect to fluctuations 
in the amplitude and the orientation of the condensate. 
We show here that the 2SC-S state proposed in this pa- 
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FIG. 1: Upper panel: Difference between the tliermodynami- 
cal potentials of two color neutral states, namely the standard 
one -with broken color symmetry- and the new color symmet- 
ric one, as functions of the baryon chemical potential /is at 
temperature T = 0. Lower panel; Eigenvalues fci, k2 — fcs 
of the curvature tensor obtained from the thermodynamical 
potential at the color-symmetric minimum. Strict positivity 
proves the stability w.r.t. arbitrary small amplitude fluctu- 
ations. Color neutrality has been imposed, but not electric 
charge neutrality. 

per fulfills these conditions. To do this, we perform a 
numerical analysis taking a phenomenologically accept- 
able set of model parameters, namely A = 653 MeV, 
GsA2 = 2.14 m, together with Go = jGs- First, 
we show that the difference between the thermodynam- 
ical potentials of the 2SC-6 and 2SC-S states is positive 



along the relevant range of baryochemical potentials, say 
350 MeV < /^b/3 < 600 MeV. The corresponding curve 
is plotted in the upper panel of Fig. 1. Second, in the 
lower panel of Fig. 1 we prove the stability of the 2SC-S 
solution by showing the strict positivity of the eigenval- 
ues fci , k2, of the curvature tensor 

'^--2aA;M-,3c_, 

derived from the thermodynamical potential in the 2SC-S 
state. Interestingly, the largest eigenvalue ki corresponds 
to fiuctuations in the "radial" direction (1,1,1), while the 
two eigenvalues denoting the curvature in the orthogonal 
"angular" directions are found to be degenerate (fc2 — 
ks). 

Note that a special situation occurs at a chemical po- 
tential /xb/3 ~ 550 MeV, where the 2SC-S state becomes 
energetically degenerate with the 2SC-6 state and simul- 
taneously the curvature in the angular directions van- 
ishes. This implies that all orientations of the conden- 
sate vector A are equivalent to each other, i.e. the cubic 
symmetry degenerates to a spherical one and the pref- 
erence of the 2SC-S state gets lost (one has in this case 
Af, = \/3As). In contrast, for all other values of /is the 
2SC-S state is preferred, owing to the penalty introduced 
for all other states which need finite color chemical po- 
tentials to achieve color neutrality. 




A = A [ MeV ] 

r g ^ J 

FIG. 2: The penalty for the thermodynamical potential, Sfl — 
qMfa _ f^jgc^^, induced by the color neutrality constraint 
in the plane of the order parameters Ab and Ar — Ag for 
jj,B ~ 1200 MeV. Since fj.3 = 0, the penalty is a function of 
fis and both vanish in the 2SC-S state where A,- = Ag = A;,. 

The behavior of this penalty 6^ in the condensate state 
space is illustrated by the contour plot shown in Fig. 2 for 
a baryochemical potential jiB = 1200 MeV chosen such 
that a maximal effect can be demonstrated, see Fig. 1. 
In this plot we consider for simplicity the plane spanned 
by the axes A^ — Ag and Ab, so that red and green 
colors are degenerate and one has /ia = 0. The penalty, 
arising from color symmetry breaking, is a function of /ig 
and vanish in the 2SC-S state for which A^ = Ag = Af, 
and /ig = 0. This figure strongly suggests that the 2SC-S 
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state is, in fact, the absolute minimum of the mean field 
thermodynamical potential. 

Another important feature of the 2SC-S state concerns 
the corresponding 12 quasiparticle modes. It turns out 
that for this state the dispersion relations and degeneracy 
factors are given by 

Et = |p|±Ms/3 [x2] , (13) 

Ei = V(b1±MB/3)2 + A2 [x4], (14) 

where A = VSA^ . This means that the system contains 
two gapless modes, just like in the case of the conven- 
tional 2SC-& state. However, in the present case the gap- 
less modes cannot be identified with the original u and d 
quarks of "blue" (unpaired) color in the original {r, 5, 6} 
color basis but arise as a combination of all three color 
states. 

In summary, we have shown in this letter that in the 
case of color neutrality the ground state of 2SC quark 
matter should be constructed in a "democratic" way, so 
that color symmetry is not broken by the choice of the 
orientation of the condensate vector in color space. For 
this state, the condition of color neutrality is fulfilled 
in a trivial way, since the penalty induced by otherwise 
necessary color chemical potentials is avoided. We have 
shown that this 2SC-S ground state is stable against fluc- 
tuations, thus the problem observed in ^3 is solved. The 



2SC-S state can serve as a starting point for considering 
hadronic correlations on the superconducting QCD vac- 
uum, where due to the entanglement of the quark color 
states in the new basis color neutral quark and diquark 
excitations arise along the radial direction besides col- 
ored excitations in the tangential plane. The conclusions 
drawn so far should not be qualitatively altered when the 
additional constraint of electrical neutrality is imposed, 
which is important for applications in compact stars. Fla- 
vor asymmetry induced by this constraint could possibly 
inhibit the formation of the 2SC state in regions of the 
neutron star matter phase diagram. However, provided 
that a phase transition to 2SC quark matter is accom- 
plished, it should be described by the 2SC-S state intro- 
duced in this paper. This issue, together with other ex- 
tensions of the present work, will be explicitly discussed 
in forthcoming publications. 
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